Vortex lattice melting, phase slips and decoupling studied by microwave dissipation in 

Bi 2 Sr2CaCu20 8+( 5 
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We have measured the microwave dissipation at 10 GHz in a Bi2Sr2CaCu2 08+i single crystal, 
between 75 and 50 K, as a function of the static magnetic field applied parallel to the c axis. We 
observe a clear-cut onset in the dissipation precisely at the melting field. In the liquid phase, for 
fields ~ 1 -2 kG, the in-plane and the c-axis resistivities deduced from microwave absorption display 
a different temperature behavior, in contrast with DC measurements. Our c-axis high frequency 
data and the DC results can be reconciled by including phase slips due to the motion of pancake 
vortices, as suggested by Koshelev. The microscopic origin of the c-axis dissipation is thus identified 
through its frequency dependence. 

PACS numbers: 47.32Cc, 74.72h, 78.70Gq 



In high temperature superconductors, thermal fluctu- 
ations play a crucial role in the behavior of the vor- 
tex system, inducing vortex-lattice phase transitions 
and new phases in the (H,T) phase diagram of the 
mixed state JO. A first-order transition correspond- 
ing to melting of the vortex lattice has been evidenced 
through a discontinuity in the entropy 0, at the same 
field H m where a step appears in the static magne- 
tization, in clean Bi 2 Sr 2 CaCu208+5 (BSCCO) ^ and 
YBa 2 Cu307_(5 (YBCO) single crystals (|. Local mea- 
surements using miniature Hall probes locate very pre- 
cisely this melting transition in BSCCO single crystals 
down to 40K S. DC measurements show that the in- 
plane resistance drops down to zero, within the experi- 
mental resolution, in YBCO untwinned crystals and in 
BSCCO at the melting transition §-|,|. The c-axis DC 
resistance exhibits a similar sharp decrease at the melt- 
ing temperature T m [5 10 1. It has been suggested that 
the latter force free (H \\ j c || c) dissipation bears some 
analogy with the voltage induced in a single Josephson 
junction by thermally activated phase slips occuring near 
T c Jll|-^8|,pt , inasmuch as BSCCO can be described as a 
stack of small area Josephson junctions JTsQ - Koshelev 
recently suggested that in-plane diffusion of pancake vor- 
tices pof is the mechanism which accounts for the similar 
temperature dependence of the ab-plane and c-axis DC 
conductivities measured above H m . In-plane dif- 

fusion induces phase slips between adjacent layers, lead- 
ing to a loss of phase coherence and thus giving rise to 
the so-called decoupling transition. It is presently still 
debated, from an experimental |3|,||,|l(|[l7 18 24| and the- 
oretical point of view |l4|,[ll|,^5j , whether decoupling oc- 
curs at the melting field or above. In this controversy, 
the c-axis dissipation mechanism plays a key role, as far 
as it may probe the phase coherence between the layers. 

In this paper, we report 10 GHz microwave dissipation 



measurements carried out in a BSCCO single crystal in 
the vicinity of the melting line, which is clearly identified 
through static magnetization measurements. The unique 
features of our experiment are: i) it locates precisely H m 
by the simultaneous onset of in-plane and out-of-plane 
microwave dissipation, ii) a quantitative estimate of both 
contributions shows that the ab plane microwave absorp- 
tion exhibits a thermally activated behavior as found in 
DC measurements, whereas the c-axis microwave absorp- 
tion, unlike the DC resistivity exhibits a different tem- 
perature behavior. This points to a frequency dependent 
resistivity associated with diffusion induced phase slips, 
as suggested in ref. pb]| . 

Both magnetization and microwave dissipation mea- 
surements were performed in the same BSCCO 
single crystal, shaped into a rectangular platelet 
(2 x 0.8 x 0.03mm 3 ), with T C =84K and a transition 
width (measured at 10 GHz) AT C ~ 0.4 K. Figure 1 
shows the result of dc magnetization measurements per- 
formed with a Quantum Design SQUID magnetometer. 
The inset displays a typical zero field cooled and field 
cooled magnetization curve at 65K as a function of mag- 
netic field. The data exhibit a "step" in the reversible 
regime, which allows to locate H m by a lower and upper 
value, as shown in the inset of fig. 1. The resulting melt- 
ing line is shown in fig. 1 between 75 K and 45K. The 
Clausius-Clapeyron equation for a magnetic first order 
transition yields the entropy change AS ~ 1 ke per vor- 
tex and per CuOi layer, consistent with other reported 
values ||. 

The experimental set-up used to measure the mi- 
crowave dissipation as a function of field (0-2 kG) and 
temperature (45-75 K) has been described elsewhere [ p6[ . 
The geometry is depicted in fig. 2. The sample is placed 
in a TE102 resonant cavity = 9.6 GHz) of an Elec- 
tron Spin Resonance (ESR) ESP 300E Briiker spectrom- 
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eter, at the maximum of the microwave magnetic field 
h\{uj) and in zero electric field. The static field H is 
applied parallel to the c axis. As hi(u>) is parallel to 
the ab plane (fig. 2), microwave currents flow within the 
ab plane and along the c axis. They generate a mag- 
netic moment m(to), and hence a macroscopic suscepti- 
bility x( w ) — Tn(u)/hi(w). Following Landau ]27| ], in 
the ohmic regime j(uj) = -Jj-^E(bj), (we have checked 
that when changing the magnitude of h\ from 5 mOe 
to 50 mOe, the response is indeed ohmic within our ex- 
perimental resolution), the energy W dissipated by the 
sample writes: 

W=^n uh 2 lX "(u)V (1) 

x"(w) is the imaginary part of the susceptibility. V 
is the volume of the sample. The measured quantity is 
<5x" = x"(H)-x"(H = 0). As explained further, x"(H = 
0) is negligible compared to %"(H), hence Sx" = x"(H). 
The absolute value (dimensionless in SI) of x" has been 
calibrated with respect to a Pbln(6%) platelet (same ge- 
ometry as the BSCCO crystal), within 20% accuracy. 

It is easy to show that as long as the depths over which 
the ab and c currents penetrate are small with respect to 
the dimensions along the same direction ("thick limit" 
approximation), W can be rewritten as: 

W « hK^abdida + R s , c d ie ) (2) 

with i=l(2) refering to position 1(2), as defined in 
fig. 2. di = 2 mm, e? 2 — 0.8 mm. R s ,k is the surface 
resistance probed by the current flowing along the k di- 
rection (fig. 2). Two measurements corresponding to two 
geometries of the sample are necessary in order to esti- 
mate R St ab and i? s ,e- Here, we simply rotate the sample 
by 90° around the c axis: e (crystal thickness) is un- 
changed, whereas d (platelet size) changes by a factor 2.5. 
It is clear from eq.2 that position 1 gives more weight to 
dissipation due to c-axis currents. 

Figure 3 shows the microwave dissipation measured as 
a function of field at 75 K and 50 K in position 1 and 2. 
At a field H* indicated by the arrow, we observe a clear 
change in the dissipation regime: the dissipation is essen- 
tially zero within our experimental resolution (~ 1.10~ 3 ) 
up to this field, then it starts to increase steadily reaching 
~ 10~ 2 at 400 G; this figure represents ~ 10 _1 x^ where 
Xn is the susceptibility in the normal state (T > T c ). 
Systematic measurements as a function of H have been 
performed from 75 K down to 50 K and display a similar 
change at fields H*(T). These fields are reported in fig.l: 
the line H*(T) coincides with the melting line determined 
by magnetization data. This establishes experimentally 
that the thermodynamic first-order transition in the vor- 
tex lattice induces a clear-cut onset in microwave dissi- 
pation. The next important result is that the dissipation 
at H > H*(T) is larger in position 1 than in position 2, 



which suggests a systematic contribution of c-axis cur- 
rents. One could think of this absorption being due to 
Josephson plasma resonance (JPR) [^9|,|o). Josephson 
plasma modes were carefully studied in the same config- 
uration (h\ parallel to the ab plane) at 35 GHz microwave 
frequency pl|: a sharp resonance in the field range 1 to 
7 kOe (depending on temperature) , is assigned to the 
longitudinal mode, whereas a much broader absorption 
is observed at much higher fields (up to 6 T). Since JPR 
frequency is expected to vary as H~ x / 2 , for 10 GHz such 
features would appear at fields roughly ten times larger, 
e.g. entirely out of our experimental field range. Two 
other remarks confirm that we can rule out JPR as the 
absorption process in our experiment. i)the low tail of 
such an absorption should already be visible in the solid 
phase, in contrast with the onset of dissipation occuring 
unambiguously at the melting field (this point has been 
checked in 5 other crystals), ii) we have swept the field 
up to IT and we observe a steadily increasing, slowly sat- 
urating signal, with no hint of an absorption peak [|30f . 
Therefore we think that we can rule out this absorption 
mechanism. We then rely on eq.2 in order to perform 
quantitative estimates. Checking the "thick limit" ap- 
proximation requires the evaluation of the surface resis- 
tance. The starting point is the surface impedance, which 
reads |§: 

Z s = R s + iX s = iujpgX (3) 

where A is the complex penetration depth. A depends 
on the London penetration depth, the normal fluid skin 
depth, and the complex effective skin depth arising from 
vortex motion J32| . 

Surface impedance measurements show that the dissi- 
pation in the London regime is X 'L ~ 10~ 4 , i.e. far below 
our experimental resolution. Therefore, we shall neglect 
the normal fluid contribution. The microwave dissipation 
due to vortices above the melting field, whether it arises 
from a flux flow mechanism |35|] or from a thermally acti- 
vated process, as we will see, [|2i|- p3| , |33"[ , can be described 
by an effective resistivity p v , (which needs not be the dc 
resistivity). We finally neglect the contribution of the 
Campbell penetration depth, since the magnetization is 
reversible above H m . Therefore the general expression 
given in ref. |32| for A reduces to: 

A 2 - A 2 - l -8l (4) 

A is the London penetration depth, and 5 V is the vortex 
induced skin depth defined as usual by 6 V — (2p v / now)? . 
The London penetration depths X a b{T) and A C (T) for 
BSCCO at T = OK are not precisely known. Quoted val- 
ues of A afc (0) vary from 1800A up to 3000 A @. A c (0) 
is still more uncertain. We have taken the value derived 
from our own measurements of the first penetration field 
of Josephson vortices, namely A c (0) = 220 fim j36|. The 
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temperature variation AX a b(T) and AA C (T) have been 
established by a cavity perturbation technique |34| . In 
the least favourable case, X C (75K) = 280 pm < d 2 — 
800 fim, and A &(75.K") = 0.6pm <C e = 20/xm. From the 
estimates of the resistivity at 2 kG described below, we 
find that <5„,a& ~ 2/xm, and 5„ iC ~ 80pm, hence again 
small with respect to e and di respectively. Therefore 
the "thick limit" approximation is appropriate. 

Combining eq.2, 3 and 4, we have computed p a b(H = 
2kOe) and p c (H = 2kOe) (we chose this field value in 
order to achieve a better accuracy; we checked that the 
results are similar for H = 1 kOe). We have reported p a b 
and p c in fig. 4 as a function of 1/T. Changing A a b(0) 
hardly changes p a b (< 10%). p c has been deduced using 
two very different values for A c (0), without changing the 
main trends: i) p a b exhibits a thermally activated be- 
havior which compares very satifactorily to dc resistivity 
measurements ^3| , ^2| j8|. This suggests that 10 GHz lies 
below the depinning frequency |35|] of pancake vortices 
in BSCCO (which to the best of our knowledge, has not 
been determined yet), ii) p c (H = 2kOc) does not dis- 
play the same temperature dependence as p a b, in contrast 
with dc resistivity data |2^||. Koshelev suggested that 
the c-axis conductivity at frequency w may be related to 
the dc conductivity through the characteristic phase slip 
frequency. The resistivity at finite frequency reads then 



p c (u) = 0) = p c (u>)- 



(5) 



The phase slip frequency value was given in feOl : uo 



ps — 

2 x 10 8 T p a bln(nr^ nax ) , where T is in K, p a b in pfl.cm, 
n is the density of mobile pancakes, and r max is the dis- 
tance over which the phase within a layer is no longer 
sensitive to pancake motion. Using our own p a {, values 
for consistency, we have computed p c (to — 0) from eq.5, 
adjusting r max so as to retrieve a similar temperature de- 
pendence for p c (uj — 0) and p a b- The result is shown in 
fig. 4: within our experimental accuracy, the temperature 
dependence of p c is now consistent with the dc data. The 
above procedure yields Infor*^) ~ 8, hence r max ~ 50, 
in units of inter vortex distance at ~ 2fcOe, if all vortices 
are mobile. Although this figure is difficult to estimate, 
the value we find seems reasonable. 

The phase slips which yield the c-axis dissipation may 
also be responsible for the decoupling process. It is 
worthwhile discussing the question of decoupling in this 
context: Daemen defines the decoupling field as the field 
where the supercurrent along the c axis vanishes p4[ . 
However, this theoretical description relies on the esti- 
mate of the phase difference between adjacent layers. 
Therefore the critical current which is calculated self- 
consistently is a local quantity. Indeed, Daemen specifies 
that at the decoupling field, vortex lines dissociate into 
a gas of pancakes. In such already suggested 



|L7j, one expects that the normal state resistance is re- 
covered. This occurs in our experiment at much higher 
field: 2000G at 75K or 3000G at 70K. Therefore there 
is an essential difference between the onset of decoupling 
(which takes place at the melting field) where phase slips 
induce dissipation but phase coherence is not lost be- 
tween adjacent layers Jl3[], and full decoupling. 



In conclusion, we have identified within a single ex- 
periment the first order melting transition, the c-axis 
resistivity and the decoupling phenomenon. Our analysis 
which accounts for the effect of frequency, is consistent 
with a recent theoretical suggestion which assigns the 
same microscopic origin to in-plane and c-axis dissipa- 
tion, namely thermally assisted phase slips. The picture 
which appears to be confirmed experimentally is that 
phase slips are induced by the thermally assisted diffu- 
sive motion of pancake vortices: the same mechanism 
might be responsible for melting, for dissipation along 
the c axis and for decoupling. 
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FIG. 1. First order melting line denned from dc magneti- 
zation measurements and from the onset of microwave dis- 
sipation. The arrows in the inset locate the field where the 
the jump of the dc magnetization starts (H m ,i ow ) and ends 
(H miU p) . 



FIG. 2. Sketch of the experimental geometry. The mi- 
crowave field hi lies parallel to the ab plane. The static field 
H is parallel to the c axis. Microwave currents (j a t and j c ) 
flow is displayed (seen from above). 



FIG. 3. Dissipation as a function of the applied field at 75K 
(a) and 50K (b) in position 1 and 2. The onset of dissipation, 
indicated by the arrow, occurs at the melting field. The ratio 
X1/X2 increases as the temperature decreases, reaching 2.4 at 
50K, which means that the dissipation due to ab currents has 
become negligible (see text). 
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FIG. 4. Resistivities p a (,(a;,H = 2kOe) (full circles) and 
p c (w,H = 2kOe) (full triangles) versus 1/T. The solid line 
shows the Arrhenius behavior of p a t- Values used for the 
London penetration depths are A a (,(0) = 2100A (see text), 
A c (0) = 220pm Q (down triangles) and 40pm Q ( up 
triangles). Open triangles display c-axis resistivity data cor- 
rected for the effect of frequency (see text). 
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